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Abstract. The Cosmic Defect theory (CD), which is presented elsewhere in this conference, introduces in the standard 
Einstein-Hilbert Lagrangian an elastic term accounting for the strain of space-time viewed as a four-dimensional physical 
continuum. In this framework the Ricci scalar acts as the kinetical term of the strain field whose potential is represented by 
the additional terms. Here we are presenting the linearised version of the theory in order to analyze its implications in the 
weak field limit. First we discuss the recovery of the Newtonian limit. We find that the typical static weak field limit imposes 
a constraint on the values of the two parameters (Lame coefficients) of the theory. Once the constraint has been implemented, 
the typical gravitational potential turns out to be Yukawa-like. The value for the Yukawa parameter is consistent with the 
constraints coming from the experimental data at the Solar system and galactic scales. We then come to the propagating 
solutions of the linearised Einstein equations in vacuo, i.e. to gravitational waves. Here, analogously with other alternative 
or extended theories of gravity, the presence of the strain field produces massive waves, where massive (in this completely 
classical context) means subluminal. Furthermore longitudinal polarization modes are allowed too, thus lending, in principle, 
a way for discriminating these waves from the plane GR ones. 
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FRAMEWORK 

The Lagrangian density we want to analyse is 

Jz? = R + -C| 1V po-e' iV e'" J , 

where the second term on the r.h.s. is a potential term, in the form of an elastic potential whose field is the strain tensor 
£/iv = \ (guv ~ i]pLv), and the Ricci scalar is computed by means of the observed metric g^ v . 

The strain is ascribed to a cosmic point-like defect, this is why we called the theory the Cosmic Defect theory [1] In 
our analogy we are interested in an isotropic medium; in this case the elastic constants take a simple form: 

Cjivpa = Ar^v^pa + Ai(7?w>*7vo- + t7mo-i?v P )- 

They now depends on two parameters only, the so-called Lame coefficients X and ji. 
The Elastic potential then translates in 

y = i[Ae 2 + 2/i£ /tv e' lv ], 

where, following the structure of the elastic coefficients, the covariant version of the e tensor is obtained by lowering 
the indices with the total metric g. 
The action integral is then 

s = jd 4 x[R + y}^. (i) 

By varying the action in eq.(l) w.r.t. the metric g^ v that is the only dynamical field we obtain the Elastic Einstein 
Equations (EEE): 

Gfiv = T^y (2) 

where 

T uv = ^ E 



8/j.v ~ £[iv 



2£jia& 



WEAK FIELD LIMIT 



In order to look for the gravitational waves of this theory we must linearize it around the Minkowski space-time, to 
which the theory reduces locally when the defect is not so strong. 
It means that the metric can be written as 

g»v = V + £ V £<Cl, 
where the perturation metric h^ v is driven by e. The strain tensor becomes 

£uv = ^ [Suv — Ifiv] = 2 £ V V ' 

where we use the coordinate invariance of the theory in order to fix the flat metric to be the Minkowski one (i.e. 
Diag(— 1,1,1,1)). The strain tensor represents the perturbation itself. Looking at the EEE, they are: 



where 



Te/xv — Ae 



e 1 

4 + 2 



(3) 
(4) 



First of all we should look at the linearized Einstein tensor, that, from 



R 



reduces to 



d r dph r a + d r d a hp y - \3h a p - d a dph - rj a pd r d s h r s + rj a pDh 



where indices are raised and lowered by means of the Minkowski metric and h is the trace of the perturbation 1 . 
The linearised Elastic energy-momentum tensor becomes 



The linearised EEE in vacuo reduce to 

[dydphl + d r d a hp y - Uh a p - d a dph -ri a pd r d 5 h\ + ri a p\3h 

A 



2\i 



hrinv 



(5) 



In order to investigate more deeply these equations we rewrite them by using their divergence and trace. The 
divergence, thanks to the contracted Bianchi identities, gives a relation between the divergence of the perturbation 
and that of its trace: 



By tracing eq.(5) one gets 



A 



2jx 



21 + — )nh-n[l + — )h = 0. 



2A 



(7) 



The trace of the perturbation becomes a dynamical field. This degree of freedom is always ghostlike, regardless the 
combination of the parameters [2, 3] (and it could be tachyonic, too): this makes the theory pathological from the 
quantum point of view. To avoid this behaviour we must constraint our parameters so that 1 + A/2ju = 0. This choice 
reduces our linearised Elastic energy-momentum tensor to the Pauli-Fierz mass term [4, 5], which actually leads to 



1 It is worthwhile here to be precise about the sign convention on Riemann and Ricci tensors. In this paper we use the Riemann tensor written 
above, where the upper index is the first one, and we contract the first and the third index in Riemann to obtain Ricci. 



different predictions from those of General Relativity, no matter how small the graviton mass is [6, 7]. This is what is 
called van Dam-Veltman-Zakharov discontinuity [8, 9], whose existence, accordingly to some authors, can be ascribed 
to the explicit breaking of the gauge invariance by the mass term so that it can be cured if the mass is generated by 
the compactification of a higher dimensional theory [10] or to the evaluation of the linearised theory around flat 
backgrounds [11]. 

The debate on massive gravity is still open and we choose the Fierz-Pauli mass term as the linearised version of the 
CD theory. 

Coming back to eq.(7), the choice A = — 2jU makes the trace of the perturbation to vanish in vacuum and now eq.(6) 
translates in what would have been the Lorentz gauge in General Relativity. This really simplifies the linearised 
Einstein tensor that reduces to — Oh^ v /2 so that the equations one finally gets are: 

(□ + M)V = (8) 
fcfv V - (9) 
h = 0. (10) 



Gravitational waves 

The set of eqs.(8, 9, 10) represent propagating massive waves and plane waves are solutions of this equations: 

V = «mv^. (11) 

Designating the propagation direction as the z axis, the wave vector is k = (to, 0, 0, ck) and from the eq.(l 1) we obtain 
the dispersion relation: 

co = ±c^/k 2 - ji 

The waves are subluminal, which is commonly referred to as being "massive" (/x < 0, consistently with the 
cosmological limit). Let us now look at which are the polarisation modes of this spacetime. From eqs.(8,9,10), we 
have 5 dynamical degrees of freedom but this does not mean that the expect the same number of polarisation modes. 
In a metric theory of gravity there can be at most six polarisation modes, as shown in [12]. The analysis can be 
performed by looking at the geodesic deviation equation, that states which is the displacement between a pair of free- 
falling particles when a gravitational wave arrives. The three-acceleration depends on the "electric" components of the 
Riemann tensor (Rq^)- It can be shown that there are six algebraically independent components of the Riemann tensor 
by using the Newman-Pensore formalism. First, one choose a complex null basis, the so-called null tetrad (k,l,m,m), 
that is related to a cartesian system {t,x,y,z} by 

k = ^=(1,0,0,1), 1 = -^=(1,0,0,-1), 

m = i=(0,M,0), m=-^(0,l,-i,0). 

We remember that we have chosen to orient the axes so that the wave travels in the +z direction, and, u being 
u = t — z/c, the k vector is proportional to Vu. Then it is possible to split the Riemann tensor into irreducible parts [ 14] : 
the Weyl tensor ( v Po, < l / i> v I'2, v P3, v P4), the traceless Ricci tensor (<£>oo><I>Oi,<I>ii. 4>i2, ^22, ^02, that are five complex 
scalars) and the Ricci scalar (A). 

When considering plane waves only some of them are different from zero and, among these, six are independent. The 
ones that are helicity (s) eigenstates under rotations about the z axis are 

X P 2 ,*22 -> s = 0, 

W 3 ,% - S=l, 
^4,^4 -» S = 2. 



We recall that these six wave amplitudes are observer dependent but there are some invariant statements that are true 
for all standard observers 2 if they are true for anyone. These statements constitute the E(2) classification of waves, 
based on the Petrov type of the Weyl tensor [13]. 

These are related to the Riemann tensor, in cartesian coordinates, as follows [12]: 

¥2 (") = -^R z o z o(u), 



*P 4 (m) = -R x0x0 (u}+RyOyo{u)+2iR x0y0 (u), 
^22 («) = - (R x q x q(u) + R }fl}! o(u)) . 



What we measure in a detection experiment is the relative acceleration of test masses, that is the six "electric" 
components of the Riemann tensor. One can express these informations in the so-called driving-force matrix 

Sij(t)=R iOj- 

In general there are eight unknowns, six polarisations and two direction cosines, but if one can establish the direction 
of the gravitational wave by other information, i.e. the k direction is known, the six elements of the Riemann tensor 
are sufficient to determine the amplitudes of the gravitational waves. 

Let us now apply this approach to our theory, where eqs.(8,9,10) must be satisfied. The second set, applied to a 
perturbation that propagates in the positive z direction, shows that the ho^ modes are proportional to the h^ z ones. The 
null trace condition, instead, allows us to express the /zoo or the h zz = (O 2 hoo/k 2 to the /?„ and h yy modes. In our case 
we still have all the six polarisation modes, as can be seen by computing the linearised Riemann tensor. 
Coming back to the driving force matrix, it can be expresse in terms of the basis polarization matrices in the z direction: 

Sy(0 = X>(z,O r «(z)y. 

r=\ 

where the amplitudes p r (z,t) are real and the index r runs over the six modes. The polarization tensor has the form 
[12]: 

e(z)h = 



*(4 



«Wy = 9 1 ' < z % = -n 1 ; (12) 



the first tensor being related to ¥2, the second and the third to the real and imaginary part of ¥3, the two next are those 
that correspond to ^4 and the last one is relative to the scalar 4>22 mode. 

These modes can play a role in discriminating among theories of gravity, and in particular they can leave a signature 
on the CMB anisotropics [15, 16] 



Static weak field limit 

When we want to reduce to Newtonian limit fields and eventual sources are taken to be static 
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2 To determine standard observers each observer sees the wave travelling in the z-dir and measures the same frequency for a monochromatic wave. 



Looking at eq.(8) we get a screened Poisson equation: 

v 2 v = 

whose solution is a Yukawa-like potential. 



(13) 



OBSERVATIONAL CONSTRAINTS 

The only parameter left in our theory is jj,. In order to quantify it we may have recourse to the fit of the type la 
supernovae luminosity which we presented in [1]. There we found a value for the bulk modulus; it was B <~ 10~ 52 
m~ 2 . From this result we have 

\H\ ~ 1(T 51 irT 2 . 

In order to compare it with the upper bounds that we find in literature it is worthwhile to rewrite our "mass" parameter 
by using the Planck constant h and the speed of light c so to get it with dimensions of a mass: 

m g = VFj- ~ 6 • l(T 66 kg. 

We know that General Relativity passes all Solar System tests so that this immediately provides an upper limit for the 
H parameter that determines the Yukawa-like fall off of eq.(13) [17]. Moreover, the absence of this effect at the galaxy 
cluster level provides the limit we were able to find [18]: 

m g <2- l(T 65 kg. 

Other limits come from the dispersion in gravitational waves since, if the graviton had a rest mass, the decay rate of 
an orbiting binary would be affected. As the decay rates of binary pulsars agree very well with GR, the errors in their 
agreements provide a limit on the graviton mass [19, 20, 21] but this limit is dramatically less restrictive than the one 
from the Yukawa potential. There are a lot of work on similar effects on the timing of a pulsar signal propagating 
in a gravitational field [22], or on the measurements of dispersion in gravitational waves using interferometers or by 
observing gravitational radiation from in-spiralling orbiting binaries [23, 24, 25, 26, 27]. 
Finally, an exhaustive review on "massive" gravitons has been done by Goldhaber and Nieto [28]. 
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